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INTRODUCTION 


When  a  slender  pointed  body  is  placed  at  a  large  angle  of  incidence,  the 
boundary  layers  on  its  sides  separate,  and  the  resulting  free  shear  layers  roil 
up  into  /ortices,  which  lie  above  the  .body  in  a  direction  between  those  of  the 
body  axis  and  the  free  stream.  These  vortices  have  an  important  influence  cu  the 
aerodynamics  of  many  combat  aircraft  and  manoeuvring  weapons. 

It  has  been  known  for  many  years  that  the  vortex  flows  of  this  kind  that 
take  place  over  sharp-edged  wings  can  be  represented  successfully  by  inviscid 
models  and  there  is  evidence  that  the  same  is  true  for  flows  over  round-edged 
wings  and  bodies,  provided  the  location  of  the  lines  on  the  surface  at  which 
separation  occurs  can  be  specified.  Since  the  noses  of  pointed  bodies  are 
approximately  conical,  knowledge  of  the  vortex  flow  over  cones  is  not  only  valu¬ 
able  in  itself  but  also  essential  for  understanding  or  calculating  the  develop¬ 
ment  of  the  vortical  flow  along  the  body.  The  consideration  of  cones  leads  tc  an 

important  simplification,  since  the  supersonic  inviscid  flow  over  a  cone  is  coni¬ 

cal,  in  the  sense  that  the  velocity  and  pressure  are  constant  along  each  ray 
through  the  apex.  Moreover,  if  the  apex  angle  is  small  enough  for  slender-bodv 
theory  to  apply,  the  lifting  component  of  the  flow  field  is  conical  at  all  Mach 
numbers . 

Use  of  the  slender-body  approximation  reduces  the  calculation  of  the  lift¬ 
ing  behaviour  to  the  solution  of  a  quasi-planar  problem  of  incompressible  flow  in 
cross-flow  planes  normal  to  the  flow  direction  or  the  body  axis.  The  vortices 
formed  by  the  rolled-up  shear  layers  can  then  be  represented  most  simply  by  the 
line  (or  point)  vortices  of  planar,  incompressible  potential  flow.  A  model  of 

this  kind  was  formulated  by  Bryson'  for  flow  past  bodies  of  revolution,  and  he 

presented  numerical  results  for  circular  cones  with  a  particular  angular  position 
of  the  separation  line  typical  of  turbulent  flow.  Additionally  he  showed  that,  as 
the  angle  of  incidence,  at  ,  is  reduced,  the  vortices  reach  the  surface  of  the 
cone  when 

-r  =  1.5  cosec  9  ,  ( 1 ) 

a  s 

where  5  is  the  semi-apex-angle  of  the  cone  and  9^  defines  the  positions  of  the 
symmetric  separation  lines  OS.  and  OS  Fig  !. 

Moore“  re-examined  the  behaviour  of  this  model  for  circular  cones.  She 

found  that,  if  9  is  smaller  than  a  critical  angle  of  about  46°,  there  are  two 
s 

solutions  for  a  range  of  values  of  a/6  smaller  than  those  given  by  (1).  For 

these  values  of  •)  ,  the  model  predicts  that  the  initial  movement  of  the  vortices 
t  s 

1 

l  _ 


4 


away  from  the  surface  of  the  cone  is  associated  with  a  decreasing  angle  of  inci¬ 
dence,  which  is  implausible  in  a  model  of  the  real  flow.  It  emerges  from  recent 

3 

calculations  t>y  Pidd  that  solutions  for  values  of  a/6  less  than  1.5  cosec  5 

s 

are  unstable  to  non-conical  asymmetric  disturbances,  so  equation  (1)  does,  in 
some  sense,  provide  a  lower  bound  to  the  values  of  a/5  for  which  solutions  of 
physical  interest  exist. 

It  i=  therefore  of  interest  to  apply  this  model  to  obtain  a  generalization 
of  (1)  to  cones  of  elliptic  cross-section .  Before  presenting  the  derivation  of 
this  generalization,  we  explain  briefly  why  we  believe  the  behaviour  of  so  crude 
a  flow  model  is  relevant  to  present-day  aerodynamics. 

Much  more  representative  models  exist  and  their  properties  can  be  deter¬ 
mined  computationally,  so  what  r61e  remains  for  simple  models?  The  more  complete 
models  that  can  be  applied  to  problems  of  vortex  flow  are: 

(a)  the  vortex-sheet  model,  which  retains  the  framework  of  slender-bcdv 
theory  but  provides  a  realistic  representation  of  the  shedding  of  circulation  and 
of  its  subsequent  convection,  in  the  limit  of  infinite  Reynolds  number; 

(b)  the  discretized  Euler  equation  model,  which  more  correctly  represents 
effects  of  chree-dimensionality  and  compressibility,  but  includes  unreal  effects 
of  numerical  diffusion,  while  excluding  the  effects  of  molecular  or  turbulent 
shear  stresses;  and 

(c)  the  discretized  Navier-Stokes  equation  model,  with  or  without  a  model 

of  turbulent  stresses,  which  includes  all  the  relevant  physical  effects  (subject 

to  the  limitations  of  the  turbulence  model),  perhaps  masked  by  non-physical 

effects  of  numerical  diffusion.  A  dramatic  justification  for  the  use  of  the 

simplest  model,  the  one  used  in  the  present  work,  is  provided  by  the  way  in  which 

it  has  contributed  to  the  understanding  of  asymmetric  vortex  flow  past  slender 

pointed  bodies  at  incidence.  The  existence  of  asymmetric  separated  flow  solutions 

of  this  model  applied  to  a  circular  cone  with  specified  symmetrically  placed 

4 

separation  lines  was  demonstrated  in  1981  .  These  solutions  were  used  m  the 

construction  of  asymmetric  vortex  sheet  solutions,  reported  in  1982^.  Solutions 
of  the  discretized  Euler  equations  for  conical  supersonic  flow,  exhibiting  pro¬ 
perties  like  those  of  the  vortex  sheet  solutions,  appeared  in  1988  .  Early 

attempts  to  find  Navier-Stokes  solutions  exhibiting  gross  asymmetry  failed,  but 
quas i-conica l ,  grossly  asymmetric  solutions  for  supersonic  speeds  appeared  in 
1989  .  All  these  solutions  are  consistent  with  a  bifurcation  of  the  symmetric 

solutions  from  the  symmetric  solution  branch  at  a  critical  value  of  the  parameter 
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a/5  ,  but  the  identification  of  the  bifurcation  locus  has  only  been  possible 

3 

using  the  present  simple  model  .  Moreover,  the  simple  model  has  allowed  the 

stability  of  the  conical  solutions  to  three-dimensional  spatial  disturbances  to 
3 

be  determined  .  With  the  vortex-sheet  model,  both  stable  and  unstable  solutions 
can  be  found,  but  so  far  it  has  not  been  possible  to  distinguish  between  then. 

The  only  solutions  that  can  be  obtained  with  the  discretized  models  are  those  for 
which  the  numerical  algorithm  is  stable,  a  condition  which  may  or  may  not  corres¬ 
pond  with  either  physical  stability  or  the  mathematical  stability  of  the  solution. 
Finally,  only  the  simple  model  gives  a  clear  indication  of  its  behaviour  near  the 
value  of  a/5  for  which  separation  first  arises,  the  particular  aspect  con¬ 
sidered  in  the  present  work. 

2  THE  MODEL 

2 . 1  Description 

We  list  the  assumptions  made  to  obtain  the  highly-simplified  model  used  in 

Refs  1  and  2  and  herein.  Firstly  we  assume  that  the  body  is  sufficiently  smooth 

and  slender,  and  the  angle  of  incidence  is  sufficiently  small,  for  slender-body 

theory  to  apply  to  the  flow  outside  the  viscous  regions  at  the  free-stream  Mach 

number  of  interest.  Secondly,  we  assume  that  the  Reynolds  number  is  high  enough 

for  the  displacement  effects  of  the  boundary  layers  and  free  shear  layers  to  be 

negligible  and  for  the  rotational  effects  of  the  free  shear  layers  which  originate 

at  separation  to  be  represented  adequately  by  vortex  sheets.  Finally  we  assume 

that  each  vortex  sheet  can  be  replaced  by  a  combination  of  a  line-vortex  and 

feeding  sheet,  as  described  below.  There  is  evidence  from  solutions  for  wing 
3  4. 

flows  ’  ,  that  this  assumption  is  reasonable  when  the  vortices  are  small  and  weak. 

According  to  slender-body  theory,  the  disturbance  velocity  potential  has  two 
components,  one  of  which  is  not  relevant  to  the  present  model  since  it  affects  only 
the  streamwise  disturbance  velocity.  The  other  component  is  obtained  as  a  solu¬ 
tion  of  the  two-dimensional  form  of  Laplace's  equation  in  the  cross-flow  plane, 
satisfying  appropriate  boundary  conditions  at  infinity  and  the  condition  of  no 
flow  through  the  body  surface.  This  component  is  itself  the  sum  of  two  parts,  one 
arising  from  the  attached  flow  around  the  body  and  the  other  representing  the 
effects  of  separation.  The  simplest  representation  of  the  effects  of  separation 
I  is  chosen,  in  which  the  axial  component  of  the  rotation  in  each  vortex  sheet  is 

concentrated  into  a  single  line-vortex,  shown  as  OV^  and  0V?  in  Figs  1  and  2. 
The  strength  of  the  vortices  is  determined  by  applying  a  generalization  of  the 
Kutta  condition,  namely,  that  the  flow  at  the  separation  line  OS  ^  or  0So  is 
directed  along  the  separation  line.  A  proper  inviscid  model  would  involve  vortex 


sheets,  rather  than  line  vortices,  and  the  pressure  would  be  continuous  across 
each  sheet.  In  our  simpler  model  we  retain  an  integral  form  or  this  condition, 
that  the  cross-flow  component  of  the  total  force  acting  on  every  elementary 
streamwise  length  of  each  'vortex  system'  vanishes.  The  vortex  system  is  the 
line  vortex,  OV  ,  and  the  cut,  or  feeding  sheet,  OSV  ,  which  joins  the  line- 
vortex  to  the  separation  line.  The  cut  is  an  essential  part  of  the  model  since 
the  strength  of  the  vortex  alters  along  its  length  in  apparent  violation  of 
Kelvin's  theorem.  We  can  regard  the  cut  as  a  feeding  vortex  sheet,  conveying 
rotation  from  the  body  to  the  line-vortex.  The  vorticity  vector  in  the  cut  lies 
in  the  cross-flow  plane,  so  the  divergence-free  character  of  the  vorticity  field 
is  maintained  without  affecting  the  cross-flow  velocity. 

2 . 2  Construction  of  velocity  field 

The  body,  an  elliptic  cone,  is  shown  in  Fig  2.  We  use  the  right-handed 
rectangular  Cartesian  axes  shown,  origin  0  at  the  cone  vertex,  Ox  along  its 
axis,  Oy  to  starboard  and  Oz  completing  the  system.  The  semi-axes  of  the 
elliptic  cross-section  are  a(x)  =  a'x  and  b(x)  =  b'x  parallel  to  the  v  and 
z  axes  respectively,  where  primes  denote  derivatives,  constant  in  the  present 
case.  The  onset  flow  has  speed  U  and  is  inclined  at  an  angle  a  to  Ox  . 

Since  a  is  small,  we  can  take  its  components  as  U  parallel  to  Ox  and  otU 
parallel  to  Oz  . 

Le  t 

Z  =  y  +  iz  (2) 

be  a  complex  variable  in  the  cross-flow  plane  x  =  const  .  Then,  if  W  is  an 
analytic  function  of  Z  ,  its  real  part  <f?{W}  is  a  solution  of  the  two-dimensional 
form  of  Laplace's  equation.  We  therefore  seek  a  disturbance  potential  in  the  form 
of  the  real  part  of  an  analytic  function  which  satisfies  appropriate  boundary 
conditions.  We  note  that  the  cross-flow  velocity  components  v  parallel  to  Oy 
and  w  parallel  to  Oz  are  given  by 


v  -  iw 


dW 

dZ  ‘ 


The  representation  of  the  onset  velocity  then  requires  that 


(3) 


W  ~  -  ictUZ 


as  Z  -*■  00  . 


(4) 


If  the  equation  of  the  body  surface  is  S(x,y,z)  =  0  ,  the  condition  or  no  flow¬ 

through  the  surface  is 

US  +  v3  +  wS  =  0  ,  (5 

A  V  2 

assuming,  as  turns  out  to  be  the  case,  that  u  ,  the  component  of  the  disturb¬ 
ance  velocity  parallel  to  Ox  ,  is  small  compared  with  U  .  The  equation  of  the 
surface  can  be  written  conveniently  as 

2  2  2  2  2  o 

S  =  b"(x)y~  +  a"(x)z"  -  a"(x)b”(x)  =  0  , 

and  so  (5)  becomes 


U(bb' y"  +  aa' z‘ 


"fab' 


,  z  ,  /  z 

b  aa  )  +  b  yv  +  a  zw  = 


(6) 


Using  o  ,  the  eccentric  angle,  as  a  parameter  on  the  elliptic  cross-section,  so 
that 

y  =  a  cos  $  ,  z  =  b  sin  i  ,  (7) 


we  find  the  body  boundary  condition  (6)  becomes 

bv  cos  o  +  aw  sin  j  =  U(ab'  sin"}  +  a'b  cos"})  .  (S) 

The  normal  direction  to  the  ellipse  in  the  cross-flow  plane  is  b  cos  :  :  a  sin  i 
so  that  v^  ,  the  component  of  the  cross-flow  velocity  normal  to  the  elliptic 
cross-section,  directed  into  the  fluid,  is 


v  =  (bv  cos  p  +  aw  sin  })/D  ; 

n 


(9) 


2.2  1  21 

D  =  (a"  sin"}  +  b~  cos"})'  >  0 


(10) 


Hence  by  (3)  and  (9) 


2  2 

v  =  U(ab  sin"}  +  a'b  cos"})/D 
n 


(1.) 


The  attached  flow  is  now  given  by  the  real  part  of  the  complex  potential 
which  satisfies  (4)  and  (11).  This  is  found  by  mapping  the  region  of  the 
Z-plane  external  to  the  ellipse  on  to  the  region  of  a  7-plane  external  to  a 
circle.  If  the  circle  is  centred  at  the  origin  and  has  radius  R  ,  a  general 
point  on  it  is 


V 


s 


Re 


10 


R(cos  ?  +  i  sin  0 ) 


12) 


so  chat  cos  o  =  i(;/'R  +  R/c)  and  sin  p  =  li(R/;  -  ;/R)  .  Hence,  if 

„  a  +  b'a-bR 
Z  "  2  R  2  ; 

the  general  point  (12)  on  the  circle  corresponds  to  the  general  point  (7)  on  the 
ellipse.  We  now  choose 

R  =  t (a  +  b)  ( i 3) 

so  that  Z  i,  at  a  large  distance,  and  obtain  the  mapping 


Z  =  ;  +  (a  -  b“)/2 


(id) 


The  normal  velocity  required  on  the  circle  in  the  ;-plane  is 


dZ  ’ 

V  =  V  — ;  , 

n_  n 1 d, , 


so  that ,  by  (9)  to  (14), 


=  U(ab'  +  a'b  +  (a'b  -  ab ' )  cos  2j)/(a  +  b) 


(15) 


Since  the  body  is  a  cone,  we  can  write 


b(x)  =  :a(x) 


(16) 


1 


where  t  is  a  constant,  the  thickness  ratio  of  the  cross-sections.  Then  (15) 
reduces  to 

v 

n  , 

'  _  raa 

~  ~  R 

This  uniform  normal  velocity  can  be  generated  by  a  source  at  the  origin  of  the 

f-plane  of  strength  2-rRv  =  2nL’xaa'  .  The  appropriate  behaviour  (4)  at  a  large 

n 

C 

distance  requires,  by  (14),  that 

W  *”  -iab'C 

and  a  complex  potential  with  this  behaviour,  generating  no  normal  velocity  on  the 
circle,  is 
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F 


-  ial' ( ;  -  R2/;) 


Henc-_  the  complex  potential  of  the  attached  flow, 
these  contributions: 

Wa  =  taa'U  log  £  ~  iaU(£ 


W  ,  is  obtained  bv  adding 
a 

•7 

-  R“/-,)  • 


To  represent  the  effects  of  separation,  we  require  the  complex  potential 
to  behave  like  that  of  a  line  vortex  of  strength  F  near  V,  ,  ie  near  Z  =  Z,  ; 
and  like  that  of  a  vortex  of  strength  ~F  near  V?  ,  i.e  near  Z  =  — Z  (  ,  where 
the  overbar  denotes  the  complex  conjugate.  The  leading  term  in  the  behaviour  of 

the  flow  at  a  large  distance  and  the  normal  velocity  on  the  body  surface  must  --t 

be  affected.  This  is  achieved  by  the  introduction  of  line-vortices  of  strength  ' 
at  ;  =  ;  and  -F  at  r  =  -r  ,  together  with  image  vortices  of  strength 

at  1  =  k }  and  F  at  r  =  — R“ / "  (  ,  where  £  j  is  related  to  Z;  through  the 

mapping  (12).  Adding  the  appropriate  complex  potentials  to  (17)  produces 

r  C  -  C,  '  +  R2/; 

+  271  i0- - ^  - 2~^  •  (!S) 

c  +  ^ —  R  /  f , 


=  'aa'U  log 


iaU  ;  - 


Positive  circulation  is  anti-clockwise  looking  upstream.  Symmetry  precludes  the 
introduction  of  a  vortex  at  the  centre  of  the  circle. 


2 . 3  Implementing  the  conditions 

We  now  wish  to  express  the  conditions  at  the  separation  lines  and  the  con¬ 
ditions  of  zero  overall  force  on  each  vortex  system,  as  described  in  subsection  2 .  1  , 
in  equation  form.  Bv  symmetry  it  is  enough  to  consider  the  starboard  half  of  the 
flow  field.  We  have  not  so  far  specified  how  the  vortex  position  and  circu¬ 

lation  F  are  to  depend  on  x  ,  nor  the  shape  of  the  separation  line  on  the  sur¬ 
face,  though  in  Fig  2  both  the  vortex  and  the  .  naration  line  are  drawn  as 
straight.  This  is  the  simplest  possibility  and  cc.L^-r.  nds  to  an  inviscid  flow 
which  is  conical  and  a  quasi-conical  boundary  layer  development.  It  agrees  quite 
closely  with  observations  of  laterally  symmetric  flow  on  conical  bodies  provided 
(a)  the  state  of  the  boundary  layer  at  separation  does  not  change  along  the  length 
and  (b)  either  the  flow  is  supersonic  or  we  are  not  too  near  the  base  of  the  cone 
if  the  flow  is  subsonic.  For  simplicity,  we  assume  the  lines  are  straight  and 
the  fLow  is  conical,  so  that 


d IF 
dx 


(  ]  91 


x 


and 


x 


!  0 


The  condition  imposed  by  Bryson  at  the  separation  line  is  that  tne  flow 
there  should  be  parallel  to  it.  A  rather  more  representative  condition  was 
investigated  by  M^ore“  and  shown  to  affect  the  solution.  However,  it  conplicat  ■< 
the  analysis,  so,  since  our  object  is  to  find  the  effect  of  non-circular  cross- 
sections,  we  retain  the  simple  form  of  the  separation  condition.  Since  the 
separation  line  is  straight,  it  can  be  specified  by  a  constant  value  of  :  : 

I  =  s  ,  v  =  a  cos  s  ,  z  =  b  sin  s  .  (20) 

The  direction  ratios  of  the  separation  line  are  then 

I  :  a'  cos  s  :  b'  sin  s 


Hence,  for  the  velocity  vector  to  be  parallel  to  this  line 

v  -  iw  =  L'a;  (cos  s  -  i-  sin  s)  ,  (2!) 


since  the  streamwise  component  of  the  disturbance  velocity  is  small  compared 
with  l'  .  The  velocity  is  found  from  (18).  Differentiating,  we  have 


*aa'  l' 


i  si'  I 


F rom  ( 3)  and  (14) 


dW 

d; 


(22) 


(23) 


dZ 

d; 


a(t  cos  s  +  i  sin  s) 


(24) 


at  tne  separation  line.  Hence  by  (21),  (23)  and  (24) 


dW  7 

'  —  =  LTaa'(t  *  i(l  -  r  )  cos  s  sin  s)  .  (25) 

d  4 

is  .  . 

Introducing  ;  =  Re  into  (22)  and  substituting  the  result  into  (25)  gives  an 
equation  which  reduces  to 

—  =  i  +  a'(l  -  ’)  sin  s  .  (26) 

1 

1  ' 


r(  *, 


to 


'  -  t; !  Re 


)  (  Re 


.  .  -  L  S  .  „  “  LS 

, )  ( Re  +  ;. )  t, Re 


The  final  condition  on  which  the  model  depends  is  LI. at  the  system  of  the 
line-vortex  and  cut  experiences  no  net  force  in  the  cross-flow  plane,  h'e  con¬ 
sider  the  transverse  forces  acting  on  an  element  of  length  of  each  component  of 
the  starboard  system,  and  represent  them  by  complex  quantities 


=  F  +  iF 

y  z 


c; 


where  F  and  F  are  the  components  parallel  to  the  y-  and  z-axes. 

y  ^ 

The  magnitude  of  the  force  per  unit  length  on  the  line-vortex  is  given  by 
the  product  of  the  density,  the  circulation,  and  th  fluid  speed  normal  to  the 
vortex.  Its  direction  is  normal  to  the  vortex  and  normal  to  the  fluid  velocity, 
in  the  sense  towards  the  side  on  which  the  speed  is  higher.  The  velocity  com¬ 
ponent  normal  to  the  vortex  arises  partly  from  the  small  inclination  of  the  vor¬ 
tex  to  the  free  stream,  giving  a  contribution  to  the  complex  velocity  of 


UZ 


CS) 


and  partly  from  the  velocity  of  the  cross-flow  at  the  position  of  the  vortex. 
The  complex  conjugate  of  this  velocity  is  given  by 


2“m  ik  '  in  lQg(z  -  V 


since  the  self-induced  velocity  of  the  vortex  is  zero.  The  only  term  in  W 
;  singular  at  Z  -  Z  ^  is  ( Z / 2 n 
in  terms  of  ;  -  we  find 


which  is  singular  at  Z  =  Zj  is  (Z/2*ri)  log(;  -  .  Using  a  Taylor  series  for 


! 


Hence  if  the  complex  cross-flow  ve'ocitv  at  the  starboard  vortex  is  V  ,  we  have 


dj.  , 

dW 

a 

- 

1  d2Z 

d;  ’1  i  ;| 

* 

dZ  | 

d: 

2^i 

T  0 

“  d.;~ 

dZ  +  -  1  +  -  2 

-  R  c  ♦  +  R_ 

(29) 


evaluated  at  :  =  \  .  The  complex  force  in  the  cross-flow  plane  per  unit  length 

of  the  starboard  vortex  is  now 


F  =  -  it Z(V  -  UZ, /x)  . 

V  ! 


(30) 
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The  jump  in  velocity  potential  across  the  starboard  cut  is  T  .  This  gives 
rise  to  a  jump  in  the  longitudinal  component  of  velocity  of  df/dx  =  7/x  ,  and  so. 
through  Bernoulli's  equation,  to  a  jump  in  pressure  of  -oLT/x  ,  since  the  cross- 
flow  component  of  velocity  is  continuous  across  the  cut.  The  complex  cross-flow 
force  per  unit  streamwise  length  of  the  cut  is  therefore  given  by 

F  =  i(Z,  -  a  cos  s  -  ib  sin  s)plT/x  ,  ( 3 1 ) 

c  1 

since  the  starboard  separation  line  is  at  2  =  a  cos  s  +  ib  sin  s  .  The  con¬ 
dition  of  zero  net  force  on  the  vortex  and  cut  is  just 

F  +  F  =  0  .  (32) 

v  c 

Introducing  (30)  and  (31)  into  (32)  leads  directly  to 


V 

U 


(2Zj  -  a  cos  s  -  ib  sin  s)/x  , 


(33) 


and  the  required  condition  follows  by  combining  this  with  (29).  The  terms  in  (29) 
follow  from  (14) 


dZ 

d; 


2  u2 

a  -  b 


d2Z 

d;2 


2  u2 

a  -  b 


and  from  ( 1 7) 


dWa 
d  C 


taa'  U  -  ictU  1  +  —r 


Hence,  after  a  little  simplification,  the  condition  is: 


/  2  2  2  \ 
(2Z|  -  a  cos  s  +  ib  sin  s)  ^4c“  -  a“  +  b  j 

4~^ 


raa  -  ia 


if 


2tU 


I  C, 


2  u2 

a  -  b 


2  2? 

lS  -  a  _+  b 

R  C  1  ( C  1  *  s  1  S1 

R2)(q  *  R2) 


(34) 
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Wich  Zj  related  to  ^  through  the  mapping  (14),  this  is  a  complex  equation 
connecting  z,  and  T  to  the  shape  and  incidence  of  the  cone.  Together  with  the 
real  equation  (26),  it  determines  the  vortex  position  and  strength. 

Before  proceeding,  we  introduce  non-dimensional  variables  and  tidy  up  the 
notation.  We  extend  the  definition  of  5  used  in  section  I  so  that  it  relates 
to  the  elliptic  cone: 


and  introduce  non-dimensional  variables  for  the  vortex  position  and  strength: 


r 

2rr6UR 


We  note  that,  by  (13)  and  (16) 


I  +  T  ’ 


so  that 


Then  the  condition  (26)  at  the  separation  line  becomes 


_ Y  (uiui  -  1  )  (oj  +  ui) _ 

.is  .  ,  is  -.  .  -is  .  .  -is  -. 
(e  -  tu)  (e  +  w)  (e  +  w)  (e  -  u)) 


a  .  „  1  -  t 
T  +  2  1 — T —  sin  s 
5  1  +  T 


ind  the  condition  (34)  of  zero  net  force  becomes 


(I  +  T) 


^[< 


(1  +  t)w  +  - - 4  -  cos  s  +  it  sin  s 


(I  ♦  t) 


a  /  I  \  .  1  —  x  uj(a)  +  u)  to 

7  W  +  -  +  iy  - 5 -  +  - - - Y -  +  - r 

'  M '  (I  +  1)0)“”  -  I  +  T  (uxj  —  I  )  (to  +1)  10  +  0)_ 


These  equations,  equivalent  to  three  real  equations,  determine  the  position  w 
and  strength  y  of  the  starboard  vortex  in  terms  of  the  axis  ratio  x  ,  the  inci¬ 
dence  parameter  a/5  ,  and  the  position,  s  ,  of  the  starboard  separation  line. 

<  can  readily  be  eliminated,  but  numerical  solution  is  required  in  general. 


i ; 


3  ANALYSIS  FOR  SMALL-SCALE  SEPARATION 

We  now  seek  a  solution  of  equations  (37)  and  (38),  assuming  that  the  vorti¬ 
ces  are  close  to  the  separation  lines,  expecting  to  find  that  the  vortices  are 

2 

weak.  We  follow  the  procedure  used  by  Moore'  for  the  circular  cone.  The  first 
step  is  to  rotate  the  axis  system  in  the  transformed  plane,  r  ,  so  that  the 
starboard  separation  line  lies  on  the  real  axis  of  the  new  plane.  Then,  if  t  is 
a  non-dimensional  complex  coordinate  of  the  starboard  vortex  in  the  new  plane,  we 
have 

t  =  e~1S£j  /R  =  uie_1S  (39) 


and  we  seek  a  solution  for  t  near  unity.  Following  Moore,  we  choose  the 
imaginary  part  of  t  as  the  small  parameter  e  ,  and  assume  the  existence  of  a 
solution  of  the  form 


t 


+  Ae  +  Be 


+  ie 


(40) 


where  the  omitted  terms  are  real  and  o(e  )  . 
We  eliminate  y  between  (37)  and  (38)  to  give 


(1  +  t)‘ 


0 


I  +  t)oj  +  - 1  -  cos  s  +  it  sin  s  (1  +  t)u>  - 


in  SJ  [( 


"I 


4t  .  a  /  \\  .fa  „  1  -  x  .  1 

- — — =•  -i  —  [w+—  +i  —  +2  - -  sms 

+  -r  )  2  3  \  <-/  y  1  +  T  J 


_ 4x 

(1  +  T) 


(e 


is 


w  is  .  -is  .  .  -is 

w)  ( e  +  ui  )(e  +oo)(e 


(uko  -  I )  (u)  +  u) 


i.r — — 

[_( 1  +  t)oo  -  1  +  t 

0)(u)  +  111)  U>  | 

+  -  2  +  - 

(wu)  -  1  )  (u)  +1)  0)  +  OJj 


(41) 


Instead  of  regarding  (41)  as  two  real  equations  for  the  real  and  imaginary  parts 
of  ui  for  given  x  and  a/i  ,  we  regard  the  equations  as  determining  the  real 
part  of  t  and  a/6  ,  for  given  Im{t}  =  e  and  x  .  So,  in  addition  to  the 
expansion  (40),  we  assume  an  expansion  of  the  form 


a 

6 


p  +  qc  +  o(e) 


(42) 
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We  note  that  e  may  be  negative,  but  Ae  in  (40)  is  positive  because  the  vortex 
lies  outside  the  circle.  Moore  obtained  a  solution  of  this  form  for  :  =  I  , 
so  it  is  reasonable  to  seek  one  for  general  values  of  x  . 

Most  of  the  factors  appearing  in  (41)  are  of  order  e^  ,  but  some  are 
smaller.  From  (39)  and  (40) 


and  so 


aj  =  e1S(l  +  (A  +  i)e  +  Be2  +  ...) 


tom  -  1  =  2Ae  +  (2B  +  A2  +  1)e2  +  . 


e1S  -  oj  =  e1Se(A  +  i  +  Be  +  ...) 


-is  -  is  .  , 

e  -  m  =  -  e  e(A-i+Be  +  ...) 


Hence,  in  order  to  obtain  the  expansion  of  (41)  to  order  e  ,  we  need  to  consider 
the  terms  0(e  ')  and  O(e^)  in  the  final  square  bracket.  We  write  these  as 


lX.  =  - - 

i  ,,  ,  2is 

(  +  i)e 


_ cos  s  +  e(A  cos  s  -  sin  s) _ _  e 

9  9  COS  S 

e(2A  +  (2B  +  A“  +  1)e)(cos  s  +  (iA  -  1)e  sm  s) 


Some  manipulation  yields 


1  +  _ t _  +  A  -  2B 

2Ae  +  2  2  .  2  ,  2 

2(t  cos  s  +  sm  s)  4A 


(x“  -  1 )  cos  s  sin  s 

xi  2  2  .2 

2(x  cos  s  +  sin  s) 


The  factor  multiplying  the  final  square  bracket  in  (41)  is  purely  imaginary,  say 
iY  ,  and  we  need  the  terms  0(e)  and  O(e^)  in  it.  Introducing  (42)  and  (43) 
we  find,  to  this  order, 


A2  +  I  (  „  1  -  t  . 

Y  =  c  - : -  cos  s  p  +  2  - -  sin  s  +  qe 

A  \ K  I  +  T  M 


(A2  +  l)2  +  2B ( 1 
2A( 1  +  A2) 
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The  first  term  on  the  right-hand  side  of  (41)  is  real  and  independent  of 
The  second  term  is,  to  the  required  order. 


i  m  ( uj  +  —  1  =  2ip  cos  s  -  2e(Ap  sin  s  +  ip  sin  s  -  iq  cos  s)  .  (47) 

O  \  U) 


The  left-hand  side  of  (41)  can  be  written  as  2L/(1  +  t)  ,  where  L  =  +  iL^ 

and,  after  some  manipulation,  we  find 


2  2  2 

L  =  2x  +  2Ae(l  -r  sin  s  +  x  (1  +  cos  s))  , 

2  2  2  2  2  2 
=  (1  -  t  )  sin  2s  +  2e(cos  s  -  2  sin  s  +  x  (sin  s  -  /  cos  s)) 


(48) 

(49) 


to  this  order  in  e  .  So,  with  (41)  written  as 


2— 2  (ViLi>  ■  — 


-  2ip  cos  s+2e(Ap  sin  s  +  ip  sin  s  -  iq  cos  s) 


(i+o 


(1  +  t) 


+  iY(Xr  +  iXj  , 


its  real  and  imaginary  parts  are 


- — -  L  =  - — j  +  2eAp  sin  s  -  YX.  , 

(1  +  x)  r  0  +  T  )  1 


(50) 


L.  =  -  2  p  cos  s  +  2e(p  sin  s  -  q  cos  s)  +  YX 


(1  +  x) 


z  1 


(51) 


Inspection  of  (50)  with  (48),  (46)  and  (45)  shows  the  0(e  )  term  vanishes 
identically.  The  O(e^)  terms  in  (51)  with  (49),  (46)  and  (44)  give,  after  some 
re-arrangement : 


(3A2  -  1)  sin  s  + 


sin  s  +  xr  1  cos  s  =  0 


(52) 


The  0(e)  terms  in  (50)  give,  making  use  of  the  previous  equations: 

- — — —  (1  +  sin's  +  t2(1  +  cos^s))  =  2Ap  sin  s  +  A~~- ^ ^  ( P  +  2  4 - —  sin  s  ]x 


(1  +  t)' 


+  i 


2  2 

(1  -  t  )  cos  s  sm  s 

2  2  .2 

x  cos  s  +  sm  s 


(53) 
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The  left-hand  side  of  (52)  consists  of  three  factors,  one  of  which  must  vanish. 

We  are  not  concert. ed  with  the  possibility  that  cos  s  =  0  ,  that  is  to  say  that 
the  separation  lines  lie  in  the  plane  of  symmetry.  If  we  suppose  the  middle  fac¬ 
tor  vanishes,  ie 

„  I  -  t 

p  =  -  2  - -  sm  s  , 

F  1  +  T 

we  find  that  (53)  reduces  to 


2  2  .2 

t  cos  s  +  sin  s  +  J  =0  , 


which  is  impossible.  Hence  we  must  choose  the  remaining  possibility: 


.2  1 


This  agrees  with  Moore's  result.  Introducing  it  into  (53)  gives 

2  2  2  2  2  2 
(1  +  x)  p  sin  s  =  2t  +  4(t  cos“s  +  sin“s) 


For  the  circular  cone,  t  =  1  ,  (55)  reduces  to 


p  =  1.5  cosec  s  , 


which  agrees  with  Moore's  result,  and  is  Bryson's  lower  bound  (1), 
The  0(e)  terms  in  (51)  are 


(1  +  t) 


2  2  2  2  2  2 
2  (c  -  2s  ♦  t  s  -  2t  c  ) 


„ ,  N  .  A2  +  1  .  o  1  -  t 

=  2 (ps  -  qc)  +  - - -  clp  +  2  [  -  -  s 


A  -  2B  -  I 


?  O 

2(t"c“  +  s-) 


(A2  +  l)2  *  2B(1  -  A2) 


where  c  and  s  have  been  written  for  cos  s  and  sin  s  .  We  see  at  once  that 
when  (54)  is  introduced,  the  term  in  q  disappears  from  (56),  so  that  q  is  not 
determined  to  this  order.  (56)  is  then  a  linear  equation  for  B  .  Manipulation 


IS 


is  helped  by  introducing 


_  2  2  2  2  2  .2 

E=tc+s  =t  cos  s  +  sin  s 


(57) 


so  that ,  from  (55) , 


ps 


(1  +  t)‘ 


2  2 

—  (t  +  2EZ)  , 


,  ^  1  ~  T 

P  +  2  s 


2E(2E  +  1) 
(1  +  t)“s 


while 


2,  2  _  2  2,  2 

(l-x)s  =  E  -  r  ,  ( 1  -  T  )c 


-  E  . 


The  solution  emerges  as 


2  /t  _  .  3A  tan  s  (2E2  3E  -  1) 
9  \E  +  E(2E  +  I) 


(58) 


When  t  =  I  ,  E  =  I  and  (58)  reduces  to 


B 


|  (2A 


tan  s  -  1) 


which  agrees  with  Moore's  result  for  the  circular  cone. 

We  find  y  ,  to  first  order  in  e  ,  from  (37)  and  the  equations  of  this 
section,  as 


Y 


2E(2E  +  1) 

(1  +  T)2 


Ac 


cot  s 


(59) 


where  E  is  given  by  (57).  As  remarked  above,  Ac  >  0  ,  since  the  vortex  is 
outside  the  body.  Hence,  since  E  >  0  ,  y  has  the  sign  of  s  ,  and,  from  (55), 
p  and  therefore  a/6  have  the  sign  of  s  .  This  means  that  the  model  represents 
weak  vortices  on  the  upper  surface  only,  and  the  starboard  vortex  then  induces 
anti-clockwise  rotation  when  viewed  from  downstream. 


At  this  stage  it  appears  there  are  solutions  for  both  positive  and  negative 
values  of  e  .  However,  there  is  a  further  condition  to  be  applied  to  the  solu¬ 
tions  found  so  far.  The  condition  which  has  been  applied  at  the  separation  line 
is  just  that  the  velocity  there  is  parallel  to  it.  This  condition  is  equally 
appropriate  to  an  attachment  line.  The  distinction  between  the  two  is  that  the 
surface  flow  converges  onto  a  separation  line  and  diverges  away  from  an  attach¬ 
ment  line  by  continuity.  This  definition  is  deceptively  simple  and  a  more  complete 
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discussion  is  planned  for  a  subsequent  report  by  the  second  and  third  authors. 
Fortunately,  when  the  vortex  lies  close  to  the  surface,  the  divergence  or  converg¬ 
ence  is  very  strong,  so  that  an  intuitive  approach  gives  the  correct  result. 

Define  v  as  the  component  of  the  velocity  in  the  cross-flow  plane  and  tangen¬ 
tial  to  the  body  cross-section,  positive  in  the  direction  of  0  increasing,  is 
anti-clockwise  as  viewed  from  downstream.  Then,  considering  the  neighbourhood  of 
an  attachment  or  separation  line,  on  which  the  velocity  vector  is  directed  along 
the  line,  we  see  that  if  v  is  increasing  with  0  the  surface  streamlines  are 
diverging,  while  if  v  is  decreasing  as  <J>  increases  the  surface  streamlines 
are  converging. 

If  we  write  the  real  part  of  the  complex  potential  W  on  the  surface  of 

the  cone  as  the  disturbance  potential  $  ,  then  d<5>/dp  will  be  related  to  the 

2  2.  . 

tangential  velocity,  and  a  negative  value  of  d  4> /d0  will  be  associated  with  a 
separation  line.  Differentiating  (18)  we  have 


dW 

d? 


raa'  U 

; 


iaU 


S  +  c  +  R 


(60) 


On  the  surface  of  the  cone  we  write  C  =  Re1*1  ,  with  £  =  Roj  and  T  =  2t5URy 

as  before,  giving: 


dW  raa' 


-2i 


U  d;  D  id 
Re 


ia(l  +  e  ~  y)  -  iiy  — r 


1  ui 

-  + 


ip  i p  -  ip  ,  -  iP  . 

e  -  oj  e  +aj  we  +  1  uie  -  1 


On  the  surface,  dC  =  iRei:>dP  ,  and  so 


J_  d_5  lj_  dWl  _  dW  d£l 

U  do  |U  dpf  "  )U  dc  dpi 


2aR  cos  0  +  >5Rv 


ip 


ip 


10  10  -  -10  -  -Ip 

e  -u)  e  +  ui  uj+e  w  —  e 


Hence,  differentiating  again. 


5  UR 


7 

Hi 

n 


„  ol  .  (  ie1  w 

2  -  sm  0  -  - 


.  lp- 
le  oi 


.  -ip  .  -1*- 

le  ui  ie  w 


,  iV  , 2  .  id  -.2  ,  -id  .2  -  -ip,.i 

(e  -  ui)  (e  +  ui)  (ui  +  e  )  (ui  -  e  ) 


(61) 
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We  now  wrice  <j>  =  s  ,  for  the  separation  (or  attachment)  line  in  (61),  and  note 
that  when  the  vortex  is  close  to  this  line  the  dominant  behaviour  is  that  of  the 
first  and  last  terms  in  the  large  bracket.  Hence,  as  e  -  0  , 


1  d2$ 

5UR  d,2 


Yl 


uue 


“IS 


we 


is 


(u  -  e~1S) 2  (w  -  e1S)2 


yi  (  22  n  9  j 

\(A  -  l)  e  (A  +  l)  £  / 


by  (43) 


4  Ay 


2  2  2 
(A  +  1)  c 


-  _  3E(2E  +  1)  cot  s 

2(1  +  t)2£ 


by  (54)  and  (59) 


By  (57)  E>0,so  $  =  s  is  a  separation  line  if  e  has  the  sign  of  s  , 
and  an  attachment  line  otherwise.  Taking  a/6  >  0  without  loss  of  generality, 
so  that  the  previous  argument  shows  s  >  0  ,  we  find  separation  corresponds  to 
c  >  0  .  This  means  the  vortex  lies  to  leeward  of  the  separation  line,  as  would 
be  expected. 


In  the  hope  of  clarifying  the  roles  of  the  solutions  with  e  positive  and 
negative,  we  use  the  formal  solution  to  express,  for  a  given  separation  line, 

■3  =  s  ,  first  the  nearby  vortex  position  corresponding  to  a  small  positive  value 
of  e  ,  and  then  the  nearby  attachment  line  associated  with  the  presence  of  the 
vortex.  Suppose  the  attachment  line  is  at  <i  =  a  .  Then,  regarding  (43)  as 
relating  the  vortex  position  w  to  the  separation  line,  we  have 


is 

e  ( i 


(A  +  i)£  +  Be" 


where,  since  e  >  0  ,  A  >  0  and  A  =  j A |  .  Now,  regarding  (43)  as  relating  the 

same  vortex  position  to  the  attachment  line,  with  a  small  parameter  e  <  0  , 
we  have 

u>  «  e^a(l  +  (A  +  i)e  +  Be2  +  ...)  , 

where,  since  £  <  0  ,  A  <  0  and  A  =  -  ) A ;  .  Hence,  equating  the  expressions 

for  m  : 
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e  (1  -  |  A  ( e  +  it  +  Be"  +  ...)  =  e  (1  +  j  A  [  e  +  ie  +  Be"  +  ...) 


Taking  logarithms  of  both  sides  and  expanding  log(!  +  x)  for  small 


x  ,  we  nave 


_  _  _  _  2  —  —  2  i  o 

ia-|A|e  +  i£  +  B£--i(-!A]£  +  ie)  +...  =  is  +  j  A  j  e  +  ie  +  Be^-jCiAle  +  it)  ‘ 


The  real  part  of  this  equation  shows  that 


e(l  +  0(e)) 


and  the  imaginary  part  shows  that 


a  =  s  +  2c ( 1  -  ] A | e)  ,  (63) 

where  j  A  j  =  1//3 

This  situation  is  sketched  n  the  plane  of  oj  in  Fig  3a  for  e  =  0. 1  . 

To  graphical  accuracy  the  term  Be^  is  negligible  for  the  separation  position 
s  =  30°  shown.  The  positions  of  the  separation  line,  attachment  line  and  vortex 
are  indicated  by  the  letters  S,  A  and  V,  and  lengths  e  and  As  are  marked. 

The  corresponding  configuration  in  the  physical  cross-flow  plane  2  is  sketched 
in  Fig  3b  for  r  =  0.5  .  The  points  S,  A  and  V  are  in  their  calculated 
positions,  but  the  separating  conical  streamline  SV  and  attaching  conical 
streamline  meeting  the  surface  at  A  are  guessed.  The  arrows  show  the  converg¬ 
ence  of  the  surface  flow  onto  S  and  divergence  from  A  . 

4  ANGLE  OF  INCIDENCE  FOR  VORTEX  FORMATION 

As  explained  in  the  Introduction,  for  circular  cones  the  value  of  a  for 

which  vortices  first  appear  at  the  surface  of  the  cone  provides  a  lower  bound  to 

the  angles  of  incidence  for  which  the  model  has  solutions,  as  long  as  the 

separation  line  is  to  leeward  of  a  critical  position,  9  —  46°  .  For  separation 

c 

positions  further  to  windward,  there  are  unstable  solutions  for  lower  angles  of 
incidence,  but  the  incidence  for  which  vortices  first  appear  at  the  surface  is 
still  of  interest. 

For  elliptic  cones  equations  (42)  and  (55)  show  that  this  angle  of  incidence 
is  given  by: 

i  ?  2  2  2 

a  _  _  _  2t“  +  4(t"  cos‘s  +  sin  s)  ,t,.\ 

-  =  p  =  - ^ -  ,  (64) 

(1  +  x)‘  sin  s 


22 


where  x  =  b/a  ,  5  =  R/x  =  Ha  +  b)/x  =  i  ( 1  +  x)a'  ,  and  che  separation 
is  y  =  a  cos  s  ,  z  =  b  sin  s  .  In  this  section  we  discuss  the  dependen 
a  on  x  and  s  .  The  elliptic  cones  include  the  circular  cone,  x  =  1 
the  flat-plate  delta  wing,  x  =  0  .  For  a  slender  circular  cone,  i  is 
semi-angle  at  the  apex,  so  that  a/5  is  a  natural  ratio  to  consider  for 
not  too  far  from  circular.  For  a  flat-plate  wing,  it  is  more  natural  to 
the  angle  of  incidence  to  the  semi-angle  of  the  wing  apex,  which  is  a' 
wing  is  slender,  so  we  also  consider  the  behaviour  of 


line 
ce  o  r 
,  and 
just  the 
cones 
relate 
if  the 


2  2  2  2  )  “ 
a  _  a  1  +  t  _  x  +  2(x  cos  s  +  sin~s 

a'  5  2  (I  +  t)  sin  s 


(65) 


Fig  4  shows  the  variation  of  a/a’  with  separation  position  for  the  circu- 

3 

lar  cone  and  the  flat  plate.  The  variations,  given  by  1.5  cosec  s  and  2  sm  s  , 
respectively,  could  hardly  be  more  different.  The  curve  for  the  circular  cone  is 
familiar,  showing  the  angle  of  incidence  needed  for  vortex  formation  increases  as 
the  separation  line  moves  to  windward.  The  curve  for  the  flat-plate  delta  is 
believed  to  be  new.  It  shows  there  are  solutions  for  all  angles  of  incidence  for 
separation  at  the  leading  edge,  s  =  0  ,  as  is  well-known,  though  for  very  small 

g 

values  of  a/a'  their  behaviour  is  unphysical  .  Since  separation  will  take 
place  at  a  sharp  leading  edge,  there  seems  little  reason  to  consider  other 

9 

values  of  s  .  However,  Wood  and  Roberts  have  shown  that  blowing  tangentially 
to  the  upper  surface  of  a  thin  delta  wing  from  a  slot  near  the  rounded  leading 
edge,  in  an  inboard  direction,  delays  the  onset  of  separation  and  vortex  forma¬ 
tion  to  larger  angles  of  incidence.  Since  the  effect  of  blowing  is  to  move  the 
separation  line  inboard,  the  behaviour  of  this  simple  model,  as  illustrated  in 
Fig  4,  is  consistent  with  their  observations. 

The  most  dramatic  difference  between  the  curves  in  Fig  4  is  at  the  small 

values  of  s  .  Equation  (64)  shows  that,  as  s  -*•  0  ,  a/5  increases  without 

limit  unless  x  is  zero,  when  a/5  -*  0  .  It  is  not  surprising  that  the  behaviour 

is  non-uniform  near  the  point  s  =  x  =  0  of  the  parameter  space,  since  the 

inclination  of  the  surface  at  the  leading  edge  changes  abruptly  through  90°  as 

x  reaches  zero.  Fig  5  shows  how  the  variation  of  a/a'  with  s  changes  as  x 

2 

increases  from  zero  to  unity.  For  the  lower  values  of  x  (x~  <  2/3)  ,  the 

curves  have  a  minimum  at  a  value  of  s  less  than  90°,  with  a  local  maximum  at 
90°.  For  larger  values  of  r  there  is  a  minimum  at  s  =  90°  . 

The  bunching  together  of  the  curves  for  s  near  90°  is  resolved  in  Fig  6, 
by  plotting  a/5  instead  of  a/a'  .  Each  curve  retains  its  shape,  but  its 
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relation  to  the  others  changes.  Also  included  in  Fig  6  are  curves  frr-  .-a lues  of 
i  greater  than  unity,  that  is  to  say,  for  cones  deeper  than  they  are  wide.  "'he 
particular  values  of  7  chosen  correspond  to  the  same  cones  as  before,  but  rolled 
through  90°.  As  x  increases  above  unity,  the  angle  of  incidence  needed  for 
vortex  formation  increases  steadily,  except  for  values  of  s  near  90°.  The  gen¬ 
eral  trend  exhibited  in  Figs  4  to  6  is  for  the  region  in  which  solutions  exist  to 
become  confined  to  the  larger  angles  of  incidence  as  the  thickness  ratio  increases 
except  for  separation  positions  well  to  leeward,  where  the  tendency  can  be 
reversed  to  a  limited  extent. 

No  attempt  is  made  to  discuss  the  results  in  terms  of  the  typical  positions 
of  separation  observed  in  laminar  or  turbulent  flow.  The  reason  for  not  doing  so 
emerges  from  a  consideration  of  Fig  7.  This  compares  vortex  sheet  and  line-vortex 
solutions  for  a  circular  cone  at  a/5  =  3  .  The  vortex  sheet  solution  is  calcu¬ 
lated  bv  the  method  of  Fiddes*^.  The  separation  line  for  the  sheet,  S  ,  is  at 

s 

s  =  0  ,  a  convenient  value  typical  of  laminar  boundary  layer  separation.  The 

sheet  extends  to  E  ,  and  is  joined  to  the  core  C  by  a  feeding  sheet  represented 
by  a  broken  line.  The  arc  labelled  ’40°  to  49°’  represents  the  vortex  positions 
in  the  line-vortex  model  for  separation  positions  40°  ^  s  <49°  .  It  is  apparent 
that,  with  s  =  45°  ,  the  Line-vortex  model  gives  almost  the  tame  vortex  position 
as  the  vortex  sheet  model  gives  with  s  =  0  .  The  non-dimensional  circulation,  '• 
of  the  line-vortex  for  s  =  45°  is  1.299,  and  the  point  B  on  the  vortex  sheet 
has  been  chosen  so  that  the  circulation  of  the  combination  of  the  vortex  core  at 
C  and  the  sheet  between  B  and  E  is  also  1.299.  The  initial  behaviour  of  the 
conical  streamline  separating  from  S  ,  at  s  =  45°  ,  is  sketched. 

This  suggests  chat  the  line-vortex  model  can  produce  solutions  resembling 
those  of  the  vortex-sheet  model  if  the  separation  position  is  displaced  to  leeward 
We  plan  to  present  more  evidence  for  this  later.  A  partial  explanation  for  this 
is  the  different  behaviour  of  the  two  models  near  the  separation  line.  For  the 
sheet  model  the  flow  leaves  the  surface  tangentially,  the  sheet  remaining  close  to 
the  body  for  some  distance  before  turning  away  from  it  and  beginning  to  wrap  arour. 
the  core.  For  the  line-vortex  model  the  separation  stream  surface  (or  conical 
streamline)  leaves  the  surface  at  right  angles.  Thus  similar  global  flow  patterns 
are  produced  by  the  two  models  if  the  separation  line  in  the  line-vortex  model  is 
displaced  to  leeward.  There  is  evidence  that  the  vortex-sheet  model  gives  vortex 
positions  in  reasonable  agreement  with  experiment  (the  occurrence  of  secondary 
separation  prevents  close  agreement)  when  observed  separation  positions  are  used, 
so  it  would  not  be  appropriate  to  use  observed  separation  positions  in  the  line- 


vortex  model. 


CONCLUSIONS 


3 

The  single  line-vortex  model  of  laterally  symmetric  separation  from  a  cir¬ 
cular  cone  at  incidence  has  been  extended  to  elliptic  cones.  Analytical  solutions 
have  been  obtained  for  the  case  in  which  the  vortices  are  weak  and  lie  close  to 
the  surface  of  the  cone.  In  particular,  the  well-known  result  of  Bryson*  for  the 
angle  of  incidence  at  which  the  vortices  first  appear  adjacent  to  specified  sep¬ 
aration  lines  on  a  circular  cone  has  been  extended  to  elliptic  cones. 

This  shows  that,  for  separation  lines  not  foo  close  to  the  leeward  generator, 
the  effect  of  reducing  the  vertical  axis  of  the  cross-section,  that  is,  of  making 
the  bod;-  mere  wing-like,  is  to  reduce  the  angle  of  incidence  at  which  the  vortices 
first  appear.  If  the  separation  lines  are  close  to  the  leeward  generator,  the 
tninner  bodies  require  large  angles  of  incidence  to  produce  vortices.  If  the 
vertical  axis  is  made  larger  than  the  horizontal  axis,  the  angle  of  incidence  fer 
the  appearance  of  vortices  increases,  unless  the  separation  lines  are  well  to 
leeward . 

In  the  limit  of  zero  thickness  ratio,  when  the  body  becomes  a  flat-plate 
delta  wing,  the  familiar  behaviour  of  the  wing  solution  is  recovered,  but  the  limit 
is  not  approached  uniformly. 


'I'M  Ae  2197 


LIST  OF  SYMBOLS 
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